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Abstract 

The eigenfunctions of the Two-Rotors Model are superpositions of states corresponding to pre¬ 
cessions of the rotors around two orthogonal axes. In the application of the model to a system of 
particles such a structure becomes a coherent entanglement of many particles. In Nuclear Physics 
such an entanglement has not been directly confirmed. I show that it is possible to come to a defi¬ 
nite conclusion about its existence by measuring the em transition probabilities for the J=3 member 
of the scissors mode rotational band and for higher excited states with intrinsic energy twice that 
of the scissors mode. The present results are relevant to single domain magnetic nanoparticles. 
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I. INTRODUCTION 


The Two-Rotors Model (TRM) describes the dynamics of two rigid bodies rotating with 
respect to each other under an attractive force around their centers of mass hxed at one 
and the same point. It was devised as a model for deformed atomic nuclei, in which case 
the rigid bodies represent the proton and neutron systems [T]. The low lying excited states 
predicted by this model were hrst observed |2] in the rare earth nucleus and then in 

all deformed atomic nuclei |3], and were called scissors modes, see Fig.l. 

By analogy similar collective excitations were predicted in several other systems |1] and as 
it is well known they have been clearly observed in Bose-Einstein condensates [S] . Moreover 
an application of the TRM to the evaluation of the magnetic susceptibility of single domain 
magnetic nanoparticles stuck in rigid matrices has given results compatible with a vast body 
of experimental data with an agreement in some cases surprisingly good [6] . 

Fig.l, however, while very suggestive, does not give a complete representation of the 
TRM states, because the TRM Hamiltonian has a double well potential and then at the 
classical level two states corresponding to the two minima. The present paper is devoted 
to the investigation of the consequences of this feature. In order to describe the problem it 
is necessary to define the model. I assume the two rotors to have axial symmetry [7]. The 
TRM Hamiltonian is then 


H = —Ll + —Lj + V 

2 Xi ^ 2X2 ^ 


( 1 ) 


where Xi, X 2 are the angular momenta, Xi,X 2 the moments of inertia of the two rotors with 
respect to the axes perpendicular to the symmetry axis and V is the potential interaction 
between them. I assume the potential to be a function of the angle between the axes of the 
rotors. Denoting this angle by 20 


U = l/(|cos(20)|). (2) 

This potential is symmetric with respect to 0 = vr/d and has two degenerate minima 
at 0 = 0,7r/2. So at the classical level the axes of the rotors will vibrate with respect 
to one another around these values of 0. At the quantum level the eigenfunctions will be 
superpositions of states describing these vibrations. To be definite let us dehne an intrinsic 
frame of axes 

t ^ 6 X Cl ^ 6 - Cl > ^ C2 + 6 _ /ox 

2sin6' ’ 2sin6' ’ 2cos0 
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FIG. 1: Scissors modes in the Two Rotors Model: the proton (p) and neutron (n) rotors process 
around the bisector of their axes. 

The eigenfunctions will be superpositions of states describing the precession of the proton 
and neutron axes around the (- and the rj- axes [1] . Such superpositions are constrained by 
the condition that independent inversions of the orientation of the proton and neutron axes 
are not observable. In general such constraint should be imposed on the absolute value of 
the wave functions|8]. In reference [1], however, they were enforced by requiring that the 
eigenfunctions, rather than their absolute value, should be invariant under these inversions, 
and we will use here this restrictive requirement. As a consequence they result to have the 
form schematically represented in Fig. 2. 

As far as I know the actual occurrence of such an entanglement has never been directly 
conhrmed in Nuclear Physics. I will discuss this feature of the TRM in detail and I will 
show how entanglement can be observed in atomic nuclei by studying the J = 3 member of 
the scissors mode rotational band and states with higher intrinsic energy. 

The determination of the eigenstates of the TRM requires the solution of the above 
mentioned constraint that until now was worked out case by case. Here I present a rather 
general and practical procedure to do it. In this way I hnd that the solution for states with 
higher intrinsic energy used in a previous investigation |9] is not unique, as I incorrectly 
assumed. 
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FIG. 2: The TRM Hamiltonian has a double well potential, the two wells corresponding to the 
precession of the rotors axes around the (- and r/-axes of the intrinsic frame. The eigenfunctions 
therefore are necessarily a superposition of the states describing such precessions. The requirement 
that they be invariant under inversion of the orientation of the neutron and proton axes, that is 
not observable, selects a definite superposition for each value of the total angular momentum. 

In Section III report the essentials of the TRM, in Section III I determine its eigenvalues 
and eigenstates, in Section IV I discuss what we can learn about entanglement in atomic 
nuclei from the existing data and new possible experiments, in Section V I compare with 
other theoretical approaches, and finally in Section VII I present our conclusions, including 
a conjecture concerning single domain magnetic nanoparticles [B]. In the Appendix I collect 
and derive some expressions of em operators. I set h = c = 1. 

II. THE TWO-ROTORS MODEL 

The TRM Hamiltonian acts on the direction cosines of the rotor axes Ci^G- These 
variables can be replaced by the Euler angles a, (3 ,7 that describe the orientation of the 
intrinsic frame plus the angle 6 . The correspondence {Ci, C 2 } = 7 , is one-to-one and 

regular for 0 < 0 < 7 r/ 2 . It is important to remember that this whole range of 6 is necessary 
for the transformation to be one-to-one m- 

Because of the axial symmetry the wave functions must satisfy the constraints 

(^1 ■ LiT = (^2 ■ T2T = 0 ( 4 ) 

In order to get analytic results these (weak) constraints on the wave functions were replaced 
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by (strong) constraints on the operators and solved in terms of 


L — Li -\- L 2 

S = L,-L2 (5) 


where 


d 

Sf = i— , S„ = — cot OLr , Sc = — tan 6*L„ . 
oO 


( 6 ) 


Using the above change of variables the TRM Hamiltonian becomes the sum of the rotational 
Hamiltonian of the two-rotors system as a whole plus an intrinsic Hamiltonian that in the 
reformulation of Ref. [lO] reads 

H=^ + H,ntr (7) 

where X = XiX 2 /(Xi -|- 12 ) and 


Hintr c2'^ 

, X 1 -X 2 


cot" 0X2 + tan 02^2 - ^ - 2 cot (20) ^ 


d 

tan 6 LcL„ — cot OL^Lc + iLc — 

oO 


+ U. 


( 8 ) 


4 X 1 X 2 

This Hamiltonian is invariant under separate inversions of the rotors axes. To dehne the 
action of such operators I must write the unit length vectors Ci)C 2 in terms of the intrinsic 
and global variables 


Cl = — sin Ot) + cos 0 C , C 2 = sin 0 p -|- cos 0 C . 


(9) 


Then the inversion operators can be represented as 

= Rc{T^)Rc{-)Re , 1-c,2 = (10) 

where i?((7r), R^(7r), R^(|) are rotation operators around the intrinsic axes. 

As I said it was assumed that such inversions should leave the wave functions invariant. 
Invariance under separate inversions is equivalent to the conditions 

Xc,Xc2^ = T (11) 

X^,T = T. (12) 

The range of 0 can be separated into two regions 

Si = s{e)s (j-G) , sii = s - 0 ) s (^ - 0 , (13) 
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where s{x) is the step function: s(a;) = l,x > 0 and zero otherwise. The 2 regions are 
obtained from each other by the reflection of 6 with respect to vr/d. It is convenient to 
introduce the notation 

«»/(«)=/(I-«)=/(«). (14) 

O O 

so that Sj= sii- With this notation V= V. 

The second term of Hintr is proportional to [ID] 9o |Xi - J2|/(4 XiX 2), where 

00 = (XC')-i (15) 


is the zero point oscillation parameter. It is therefore negligible for atomic nuclei (for 
which |Xi — X2|/(4 XiX 2) << 1 and 0o ~ 0.1) but not for free nanoparticles (for which 
|Xi-X2|/(4I, X 2 ) ~ 1 and 00 < !)• I think, however, that the importance of the second term 
of Hi for different moments of inertia is due to the fact that the intrinsic frame I choused is 
not a principal frame, namely a frame in which the tensor of the moment of inertia of the 
two-rotors system is diagonal. I conjecture that in a principal system the second term will 
be small also for nanoparticles. 

Neglecting the second term the TRM Hamiltonian becomes then invariant also under the 
transformation 

R = ^-) Rg . (16) 

Next I eliminate the linear derivative in the first term of Hi^tr by the transformation 


getting 


(H$)(0) = 


-^2 sin(20) 


$'( 0 ). 


= - 


+ cot^ 0 -I- tan^ 0L^ 


X_ (2 + coe(2»)) 


+ v(e). 


(17) 


( 18 ) 


At last I assume that the wave functions have such a fast falloff (which is completely justihed 
in the case of nuclei) that I can perform the harmonic approximation for the potential and 
the circular functions 


V 6^ , S'^ = iVe, = , in region I 

, S^ = -iVo, = S'^ = 0, in region II 
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(19) 

( 20 ) 






where 


I then write accordingly 


d 1 

Vfl =- 

^ de 29 


H'intr ~ HjSj + HjjSjl 


(21) 


( 22 ) 


where 


H, = -w 


Hii — -u 




d X X 


TT 






TT 


° W„ 


(23) 


with 


9 o 9 

TJ 




(24) 


(25) 


The harmonic approximation makes more evident that (18) is a double well Hamiltonian. 


implying that in stationary states the rotor axes oscillate simultaneously around the (- and 


p-axes. 


The eigenfunctions and eigenvalues of Hj are 110 


n\ 


^Kn\^) — 


^nK — 


n + K)\9o 
oj{2n + K + 1) 


X 


K+^ tK 


i LS e-i’ 


(26) 

(27) 


where are Laguerre polynomials and the wave functions ifKn are normalized according 
to 

dx {ipKn{x))‘^ = ^ • (28) 

Because in the harmonic approximation 9 plays the role a a radius, I call n the radial 
quantum number. 

In general the eigenstates occur in doublets, whose energy splitting can be estimated with 
the WKB approximation 

re{E) 

6 E^Eexp / (-|p(0)|) (29) 

J-e(E) 
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where 0{E) is the angle of inversion of the classical trajectory of energy E and p{6) its 
conjngate momentnm, \p\ = ^J\2I{E — V)\ ~ \ sm.6\/6Q. Becanse 0{E) k, 9q for the states 
of interest 

5E ^ Eex.p(^‘^ . (30) 

For atomic nnclei in the rare earth region 9‘^ ~ 0.01 and snch energy splitting is to all effects 
negligible, bnt the sitnation is different for nanoparticles. 


III. EIGENSTATES 


I write the eigenfnnctions in the form 

where 


K>0 


-pi — 
MK — 


21 + 1 


{'^MK + ( ^Y'^M-k) ■ 


(31) 


(32) 


16(1 + 6 ko)'^‘^ 

I, M, K are the nnclens angnlar momentnm and its component on the 2 ;-axis of the laboratory 
frame and the C-axis of the intrinsic frame, and m an additional qnantnm nnmber to be 
specihed in the seqnel. Becanse all the states I will consider have positive parity I will omit 
the parity qnantnm nnmber. The combination of rotational matrices in the E is reqnired 
by the condition (0- It remains to impose the condition ( [I^ . 

The eigenstates are normalized according to 


/* 2 tt /*7r /‘27r 

/ da d/dsin/d / dy / dd 
'0 Jo Jo Jo 


= 1 . 


The eigenstates of the Hamiltonian in region I are 

'^LM,K,n = J^MKi(^,^E)+KniO) . 

For each snch eigenstate there is in region II the degenerate eigenstate 

,(//) 


^/,Ar,X,n ~ Gm,K ^K,n 


where 




M,K 


(33) 


(34) 


(35) 


( 36 ) 


The constraint (12) determines the their amplitndes in the total eigenfnnction. 







When I express the Ql^ ^ in terms of the total eigenfunctions take the standard 


form (31). Notice that in region I the eigenstates have a unigue component of while in 
region II they have all the components of Iappearing in guantum number m is 

the component of the total angular momentum on the (-axis in region I. Even if each of the 
rotors has axial symmetry, the two-rotors system does not have it, so that the component 
of angular momentum along any intrinsic axis is not conserved, resulting in a superposition 
of states with different i^-quantum number. 


In order to impose the constraint (12) I must determine the action of i?,j(7r) and R^{ti/2) 
on the IFfj ^ and the Qfj For any component of /*,, k = (,7],( 


so that 


exp{ialk) = i— sin(4a) -f cos(Jfca) 

h 


(37) 


exp(i7r4)V’4 = 

exp(i7r/2 4)^/’4 = sin(47r/2)cos(47r/2) 

L J-k 




(38) 


Notice that the transformations in the last equation are simpler for 4 even. 

In order to find the action of R^{7i/2) on the Rm k the Gm x ^ express these functions 
in terms of the eigenstates of 


j2 _ t^2p^I 


(39) 


Because, as noted above, such an action is simpler for even values of K it is convenient to 
express all the Rfj ^ the Gfj ^ for K even and odd, in terms of the ^ with even K. 


IV. THE SCISSORS MODE ROTATIONAL BAND 

For the discussion of entanglement it is necessary to separate the contributions coming 
from regions I and II. To this end I introduce the parameters r/,r// that in the TRM take 
the values 

r/ = r// = 1, in the TRM . (40) 

A general feature is that the intraband magnetic transition amplitudes vanish, because they 
are proportional to 

< </2i,o|Ve|</9i,o >=< >= 0. (41) 
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A. The band head 


The band head, the scissors mode, is a pure K = 1 state. Its wave function and transition 
amplitude are well known [T] but are reported for the sake of completeness 




where 


The transition amplitude to the ground state is 


(42) 


(43) 


< ‘I'iMi,o|244(Ml;/i)|'I'o,o,o,o >— ^2W(Ml)Co,^^(r; + r;/) (44) 


where the expression of 


is riderived in the Appendix. 


Af(Ml) = \ 

V dvr 2 m 


(45) 


B. The J=2 member of the band 

The J = 2 member of the band is also a pure A' = 1 state, and its wave function and 
transition amplitude are also well known jl] but are reported for the sake of completeness 

where 

*^1,1,0 = 9^1,0+ • (47) 

Its transition amplitude to the ground state is 

< 4/2M,i,o|Ad(772;0,0,0 >= —^e(52o^^o C'o(^^ (r/+ r//) (48) 

where Q 20 is the quadrupole moment in the intrinsic frame. 


C. The J=3 member 

The wave function of the J = 3 member is determined in the present paper. It can be 
written 

4'3Vfi,o = c-Affi (fifl + s Q^ii 9 ^ 1,0 5 = 1 (49) 
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where 


~ 


ms) 


(50) 


The eigenfunctions of /| with eigenvalues 0,4 respectively are 


^Mo ~ 

“ 2^2 


M3; 


(75 + Vs J-) 


ms) 


(51) 


Expressing and in terms of K\jq and lC\j2 and imposing the constraint (12) I get 


*^1,1,0 - ‘dpo + ^ ‘dl,( 

^3 ^ ° 

*^1,3,0 - —j— ‘^1,0 • 


Written in the standard form (31) 


T -r- l / 1 ° \ a/T 5 _o o 

'h3Ml,0 = -^1 (V^io + ^ ‘dio j H— 3 <^10 


(52) 


(53) 


One can see how the intrinsic structure of the two-rotors system changes in the band with 
the angular momentum, with a strong departure from a rigid rotor. 

The nonvanishing electromagnetic transition amplitudes are 

< ^3M1,o|-^(-^ 2;/i)|\h2M'l,0 > = C'2M'2/i < ‘dlol'dlO > 

< ^3Mi,o|A 4(E2;/i)|\hiM'i,o > = eQ 20 Y^C'im'2a( <‘/^io7io > ^0.63r/— 0.5r//j 

< 4/37^1^01-^(1143;/x)|\['o, 0 , 0,0 > = < (75io|V|(poo > (^rj — 

(54) 


where the expression of 




i-h^ 

3 Vi?3 


-d 


(55) 


is derived in the Appendix. 


V. OVERTONES 

In previous papers [9III]! studied the states of intrinsic energy 2 Escissors, called first over¬ 
tones because of the harmonic approximation. I know that in general in Nuclear Physics 
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we can trust collective models at most for the lowest excitation. Nevertheless I considered 
worth while investigating the hrst overtones for two reasons. First their excitation energy 
falls below the threshold for neutron emission and therefore their width is of purely electro¬ 
magnetic nature, which might make their existence plausible, in spite of the fragmentation of 
the scissors mode. Second, their electric quadrupole transition amplitude is of zero orderjH] 
in 6 * 0 , and therefore much greater than that of the J = 2 member of the scissors rotational 
band that is of order Oq. 


I reconsider now these states by using the present method of solving the constraint (12). 


The state / = 0 = m = 0,n = l cannot be excited by electromagnetic radiation, 
and for this reason it was called the elusive overtone m- The same is true for the state 
/ = 1 , m = 0 , n = 1 . 

The states J = 2,m = 0 ,? 7 , = l and I = 2,m = 2,n = 0 are degenerate and their wave 
functions in regions I and II are 


Si d'2,M,2,0 = J^M2 ^“2,0 


(56) 


sii't 


II 

2M2,0 


— ^M2 ^ 2,0 




2M0 


,1 ~ ^MO ‘^O,! 


(57) 


where 




M 2 ) 


Gm2 — ^M2) 


(58) 


The eigenfunctions of J| with eigenvalues 0,4 respectively are 


r2 


r2 

^M2 


~ 2(^mo i 


M 2 ) 




M2) 


(59) 


Expressing the Jy 


2 

MK 


and 


MK 


m 


terms of /C|^q and ICj ^2 imposing the constraint ( 12 ) 


12 




I find 


2 1 ° 

*^ 0 , 0,1 = ^ 0,1 - 2 ‘^ O ,! 

2 ° 

*^0,2,1 =- Y 

^2 o 

*^ 2 , 0,0 -- Y 

*^^ 2,0 = + 2 ^^ 2,0 • ( 60 ) 


The state T 2 mo,i can be regarded as a member of the rotational band over the elusive over¬ 
tone To, 0 , 0 , 1 - Ifs quadrupole transition amplitude to the ground statae vanishes. Because, 
as the bandhead, it cannot be excited from the ground state, I will not discuss it any further 
(even though if reached from above, it could decay to the scissors mode). 

The nonvanishing electromagnetic transition amplitudes of the state T 2 M 2 ,o are 


< T2Ar2,o|-^(-E'2;/i)|To,o,o,o > 

< T2M2,o|A^(Tfl;/i)|Ti,M'i,o > 

The transition strengths are 






(r/ rn) ■ 


B{E2) tovertone= ^scissors 

B(M1; overtone —)■ scissors) = - B{M1) tscissors • 


(61) 


(62) 


In the quoted investigation of overtones [9] I did not hnd the present (60), most general 


solution of the constraint (12), but the particular solution 


1 

71 


(T 2 Ar 2,0 


T2M0,i) 


(63) 


which was incorrectly assumed to be unique, and the em transition amplitudes were evalu¬ 
ated accordingly. I notice that with the TRM values of the parameters r/,r//, the electric 
quadrupole transition strength of the overtone T 2 M 2 ,o is a factor 2 larger than that of the 
above state while the magnetic dipole transition strength is a factor 4/7 smaller. 


VI. ENTANGLEMENT 

The TRM gives distinctive predictions that should enable us to reach a dehnite conclusion 
concerning the existence of entanglement in atomic nuclei. 
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To be definite I compare the predictions of the TRM with those of a Reference Model 
that does not have entanglement. This Reference Model is what is often regarded to be the 
TRM as derived from microscopic models [T2HI1] . It is the intrinsic Hamiltonian Hj with the 
understanding that it acts on intrinsic wave functions defined and normalized in the whole 
range 0 < 9 < tt 12. The Reference Model has the same eigenvalues as the TRM. Unlike the 
TRM it has axial symmetry and obviously describes a precession around the ('-axis only. Its 
eigenfunctions are pure R'-states and can be obtained from the eigenfunctions of the TRM 

O 

setting ^Kn= 0 and < '^Knl'^Kn >= 1- The transition amplitudes can be obtained from the 
expressions relative to the TRM by setting 

rj = 2 , rjj = 0 , in the Reference Model. (64) 

I discuss the entanglement separately for the different states. 

A. Entanglement in the scissors rotational band 

1. Entanglement in the J=l,2 members of the band 

The em transition amplitudes of these states are the same in the TRM and in the Ref¬ 
erence Model, because they are proportional to rj + rn, a quantity that takes the same 
value in both models. Therefore we cannot learn anything about entanglement from their 
comparison with experiment. 

2. Entanglement in the J=3 member of the band 

Let us denote by RM,TRM the transition amplitudes for the Reference Model and the 
TRM respectively. Then 

< ^3mi,o|-^(-E' 2;/i)|'I'2M'i,o = —8 < T3 mi,o|-^(-E'2; 

< 4/3mi,o|-^(-E'2 ;= 9.7 < \I'3mi,o|7W(U2; 

< 4/3Mi,o|7W(M3; 0 , 0,0 ^ < ^3Mi,o| 7W(M3;/i)|\ko,o,o,o • (65) 

One can see that the amplitudes for decay of the J = 3,m = l,?7, = 0 state to the lower 
members of the band are depressed by large factors in the TRM with respect to the Reference 
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Model. Eqs.(54) show that this is due to destructive interference between the contributions 


from regions I and II. The difference in strengths is so large that if this member of the band 
can be observed one should be able to reach a dehnite conclusion about entanglement. 


B. Entanglement in the first overtones 

The electric quadrupole amplitude for decay of the overtone T 2 m, 2 ,o to the ground state 
vanishes in the absence of entanglement. The relation between the magnetic dipole transition 
amplitudes in the Reference Model and the TRM is 


< '^2M2,0 |7W(Ml;/r)|To, 0 , 0,0 1.7 < T2M2,o|-M(M1; /i)|To, 0 , 0,0 > 


TRM 


( 66 ) 


Observation of the magnetic transition in the absence of the electric decay would give strong 
support to the absence of entanglement. Obviously on the contrary, observation of both 


transitions with the strengths (62) would be evidence in favor of it. 

In a recent experiment the deformed nucleus ^®®Gd, where the scissors mode has been 
discovered initially [2] , has been studied by a high resolution nuclear resonance fluorescence 
experiment at the S-DALINAC up to 7 MeV of excitation energy. ” A single candidate with 
the following characteristics a) a ground state decay indicating a quadrupole radiation, and 
b) simultaneously a signihcant branch to the main fragment of the scissors mode at 3 MeV 
has not been found above the detection limit ” [15] . 

For an assessment of the realization in nature of the hrst overtone and its entanglement it 
is crucial to put the above Endings in relation with the present estimate of its decay strength 
to the scissors mode. Indeed such a strength is not so large and in the comparison with 
experiment it should be reduced by a factor equal to the percentage of the total strength 
carried by the main fragment of the scissors mode. 


VII. OTHER THEORETICAL APPROACHES 

There is a copious literature on the scissors modes, in which however entanglement never 
appears explicitly. Therefore it is sufficient for me to examine schematically how could one 
investigate entanglement in the different approaches. For this purposes I can schematically 
divide them into two categories. 
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In the first one, following different procedures, one derives a collective Hamiltonian that 
has an eigenstate with the quantum numbers of the scissors mode and approximately the 
same excitation energy. There can be however some important subtleties that I illustrate by 
two examples. One is provided by the Interacting Boson Model [I6] . It has been shown na 
that in the coherent states approximation, for small vibrations of the rotor axes around 
the C-axis, it reproduces the intrinsic part Hj of the TRM Hamiltonian. I think that the 
vibrations around the //-axis are also present in the IBA Hamiltonian, and that in the co¬ 
herent state approximation they should provide the Hamiltonian Hjj, but this remains to 
be verihed. I must notice, however, that in calculations done with the IBM one does not 
use the coherent state approximation, but rather other approximations assuming explicit 
symmetries of the wave functions. In a comparison with the TRM one has to check whether 
and how the invariance under inversion of the orientation of the rotors axes has been implic¬ 
itly implemented, and whether the assumed symmetries imply, for instance, axial symmetry, 
that would eliminate the entanglement altogether. 

Another relevant example is the recent analytical approach to rotational states [H], in 
which the TRM Hamiltonian has been derived in the form Q. This paper is especially 
interesting in our context, because in the derivation of the collective Hamiltonian, as far as 
I understand, entanglement has not been enforced explicitly, and then also the Hamiltonian 
of the Reference Model should be a possible outcome. A clarihcation of this point is of 
the highest consequence for a strict connection between a many-body Hamiltonian and the 
TRM. 

In conclusion one must be sure of which conditions concerning invariance under inversion 
of the rotor axes are explicitly or implicitly set on the wave functions in the course of the 
derivation. 

The second category includes model or microscopic calculations in which a collective state 
appears that can be interpreted as the scissors mode. The RPA for instance, reproduces at 
a semiquantitative level the eigenvalues and the em strengths of the TRM for scissors modes 
czi- A recent approach, the Wigner Function Moments method ra. also belongs to this 
class. 

In all the works belonging to this category, however, to our knowledge the resulting 
collective modes have not been analyzed in relation to the entanglement. 

All the theoretical approaches of which I am aware are restricted to the lowest scissors 
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excitation. This is justified by the fact that in general collective models in nuclear physics 
can be trusted at most for the first excited state. I notice, however, that this does not need to 
be an absolute rule, and indeed it is not true for all systems. For instance in the evaluation 
of the magnetic susceptibility of single domain magnetic nanoparticles using the TRM all the 
excited states appear and contribute [6] . The important point is whether the rotors actually 
behave as rigid bodies at the energy of the collective state of interest, namely whether 
the coupling between intrinsic and collective degrees of freedom is or is not important. 
A general criterion can be found in [TH]. But for higher states this point can be more 
efficiently investigated in a constructive way, introducing in a many-body Hamiltonian a 
number of collective variables with an equal number of constraints in order not to change the 
effective number of degrees of freedom. In a variant of such a method one can avoid explicit 
constraints that make the calculations akword by modifying the microscopic Hamiltonian 
in such a way as to push the spurious excitations associated with the redundant variables 
out the part of the spectrum one is interested in. Such a method has been used long ago to 
enforce translational invariance [12] in shell model calculations and exploited to introduce 
collective rotations IZDl. The latter application might be extended to the physics of the 
TRM by introducing the collective variable 6 in addition to the Euler angles. 


VIII. CONCLUSION 

The wave functions of the TRM have a peculiar entanglement. In applications of the 
model to many-body systems this becomes a coherent entanglement of many particles of 
which I do not know other examples. 

In nuclear physics with the present data there is no evidence in favor or against it, and 
the only check I can envisage is to compare the mass density distribution of the states in 
which the scissors mode is fragmented with that predicted by the TRM. 

I have shown, however, that significant pieces of information can be obtained from the 
study of higher excited states. I hope that a definitive assessment concerning the first 
overtone will come soon [T3] . The other crucial investigation concerns the J = 3 member of 
the scissors rotational band. If such a state is realized in nature and can be observed one 
has enough distinctive predictions to identify it. 

It is interesting to consider the application of the TRM to single domain magnetic 
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nanoparticles. These objects consist of a magnetic structure, called macrospin, that ro¬ 
tates with respect to a nonmagnetic lattice. They have been represented as a couple of rigid 
rotors, one associated with the nonmagnetic lattice, and the other one, with a spin attached, 
with the macrospin | 6 ]. The macrospin has usually two stable orientations antiparallel to 
each other, separated by an energy barrier. At hnite temperature there is a hnite probabil¬ 
ity for the magnetization to flip and reverse its orientation. The double well potential, at 
variance with the case of atomic nuclei in which it might appear an artifact, is in this case 
at the basis of the dynamics. There is a strong, even though indirect evidence of the validity 
of the TRM for nanoparticles stuck in rigid matrices [B]. I think that a direct check of the 
entanglement predicted by the TRM is possible by measuring the magnetic susceptibility of 
free nanoparticles at temperatures of the order of 1 K. 


Appendix A: Electromagnetic operators 


The magnetic multipole operator in the intrinsic frame is 


Al'(M/,/i) = — 


dfSi 


d 




(Al) 


TTT'C Vfiucleus J 

where Vnucieus is the nuclear volume and the operators in the intrinsic frame are given in 


Eqs. (20). I found that the terms S'^, S'^ do not contribute to the transitions of the states I 
consider, and I will ignore them. Working out the above equation I then get the expression 
of the magnetic dipole (already well known) and octupole operators in the laboratory frame 




Vo sji) 


(A 2 ) 


where 




i--f^ 

3 Vi ?3 




2 mc 


(A3) 


R 3 , Ri are the lengths of the semiaxes of the ellipsoids. In the evaluation of transition 
amplitudes I will need the matrix elements 


< ^2o\'^e\^io >— 




(A4) 
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The electric quadrupole operator in the laboratory frame was evaluated in [T1I9] 


= eQ 


20 




Si - -sii\ + 7: + T>1_2) sii 


^20 {Osi+ 6 sii) 


(AS) 


where eQ 2 o is the electric quadrupole moment in the intrinsic frame. Notice that the first 
line is of zero order in 9 while the second line is of order 9. In the evaluation of transition 
amplitudes I will need the matrix element 

1 


< > — 


2^2 


(A6) 
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